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TECHNICAL MEMORANDUM 


AUTOMATIC MATHEMATICAL MODELING FOR REAL TIME 
SIMULATION PROGRAM 

I. SUMMARY 


This report describes a methodology for automatic mathematical modeling and generating 
simulation models. The models will be verified by running in a test environment using standard 
profiles with the results compared against known results. The major objective is to create a user- 
friendly environment for engineers to design, maintain, and verify their models; to automatically 
convert the mathematical models into conventional code for computation; and, finally, to document 
the model automatically. 

The traditional method for simulation modeling is to define, derive, and organize the 
equations and then to develop the program for manual computation. If any modification or correc- 
tion needs to be done in the design or mathematical equations, it will create a tremendous amount 
of work for the rest of the process. Generally, any model will require several modifications before 
it matches a real system. Historically, the modifications have been time consuming and a fertile 
source of errors. 

Some of the available symbolic mathematics tools will help to derive the equation symboli- 
cally, and automatically produce the final equation and program code. The use of the LISP langu- 
age can build the user-friendly interface and generate the knowledge base for the symbolics mathe- 
matics tool to build equations and program codes. After the equations and codes have been stored 
in the files, the LISP program can then combine and organize them and create a complete 
compiled source code for any conventional language requirement. A demonstration program was 
designed for modeling the Real-Time Simulation system. It is written in LISP and MACSYMA and 
runs on a Symbolics 3670 LISP Machine. It contains an initial set of elements for the Real-Time 
Simulation system and a questionnaire that allows the engineer to answer a set of questions to 
specify a particular model. The system is then able to automatically generate the model and 
FORTRAN code. The future goal which is under development is to download the FORTRAN code 
to the VAX/VMS system for conventional computation. The mathematical model will be verified in 
the test environment and the solution compared with the real data profile. Control, display and 
ancillary systems will be developed which will allow the execution of a predefined profile (from a 
file) or interactive modification of the variables. Variables can be displayed in a text format or as a 
plot from the control and display system. An execution control module will be available to allow 
time and sequence control of the model and peripheral models. A recording module will be avail- 
able for execution in series with the math model and will record (optional) interface variables trans- 
parently. This module will use packet definition variables and circular queues to determine what 
and when to record. 

The use of Artificial Intelligence (AI) techniques has shown that the process of simulation 
modeling can be simplified. 
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II. INTRODUCTION 


Mathematical modeling for real-time simulation programs is a very complicated process 
which includes analysis, design, and the generation of complex equations and programs. Generally, 
the model will require several modifications before it will match a real system. Historically, the 
modifications have been time consuming and a fertile source of errors. 

The use of A1 techniques has shown that this process can be simplified. Some of the AI 
software tools will allow us to create a user-friendly environment for engineers to design, maintain, 
and verify their model and also to automatically convert the mathematical model into a conven- 
tional programming language for execution. 


III. PURPOSE 


The major objective is to develop and create a very comfortable environment for engineers 
to design and maintain their model. The automatic Mathematical Modeling generates knowledge 
base, mathematical equations, and conventional program codes. It helps to simplify the process of 
the modeling, and cuts down the development time and errors. 

The automatic generation of a math model in simulations will add confidence in the simula- 
tion. Using interfaces which are well defined and validated will reduce programming and debug 
time and will allow concentration on logic and equations. 

This technique can be applied to many different science and engineering projects. 


IV. INTELLIGENT INTERFACES 


The intelligent interfaces include: 

1) The user interface 

2) The interface between knowledge base and automatic model generation 

3) The interface between modules comprising the simulation. 

Automatic modeling requires a very friendly work environment to collect the necessary 
information and to generate the knowledge base. Through the knowledge base information, the 
system automatically creates the mathematical equations and generic program codes for equations. 
Next, all the equation codes are generated and linked by another AI program to combine and 
organize the codes together and to build a completely compiled program. 
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An intermodule interface will become part of the automatic model generation process. Vari- 
ables will be identified and an external interface will provide the variable and its type to external 
modules. These external modules will record, display, or modify the variable. 


V. AUTOMATIC MODELING 


The traditional method for simulation modeling is to manually define, derive, and organize 
the equations and then to develop the program for computation. If any modifications or corrections 
need to be done in the design or mathematical equations, it will create a tremendous amount of 
work for the rest of the process. Generally the model will require several modifications before it 
matches a real system. Historically, the modifications have been time consuming and a fertile 
source of errors. 

Some of the available symbolic mathematics tools will help us to derive the equation 
symbolically and to automatically produce the final equation and program code. The use of LISP 
language can build the user-friendly interface and generate the knowledge base for the symbolics 
mathetmatics tool to build equations and program codes. After the equations and codes have been 
stored in the files, the LISP program can then combine and organize them and create a complete 
compiled source code for any conventional language you require. 

The use of AI techniques has shown that the process of simulation modeling can be simpli- 
fied. 


Control, display, and ancillary systems will be developed which will allow the execution of 
a predefined profile (from a file) or interactive modification of the variables. Variables can be dis- 
played in a text format or as a plot from the control and display system. An execution control 
module will be available to allow time and sequence control of the model and peripheral models. 

A recording module will be available for execution in series with the math model and will record 
(or not record) interface variables transparently. This module will use packet definition variables 
and circular queues to determine what and when to record. 


VI. SYSTEM OVERVIEW 
A. Flow Diagram 
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B. Basic Software and Hardware Requirements 


The Basic software tools used were “MACSYMA” and “LISP.” The hardware was 
“Symbolics 3670” and “VAX.” We used a Symbolics for the equation and conventional code 
generation and then downloaded to the VAX system for computation. 

There are a number of tools available for Symbolic Mathematical equation and conventional 
program code generation. Some examples are as follows: 


Symbolics 
Math Tool 

Memory 

Requirements 

Available Computers 

MACSYMA 

2 Megabytes 

DEC 10 or 20 
VAX 

Honeywell 6000 Series 

Symbolics 

LM2 or 3600 

SMP 

2 Megabytes 

VAX 730, 750, or 780 
Using UNIX operating system 

REDUCE 

350 kilobytes 

IBM-360 or 370 
DEC 10 or 20 
VAX 

UNI VAC 1100 
CDC CYBER 
CRAY-1 

BURROUGHS 6700 

APOLLO DOMAIN AND OTHERS 

MuMATH 

512 kilobytes 

Personal Computer 

MAPLE 

350 kilobytes 

VAX Using UNIX 4.3 

ALTRAN 

270 kilobytes 

1966 Standard FORTRAN 

FORMAC 

150 kilobytes 

IBM-360 or compatible machine 

SAC-2 

120 kilobytes 

1966 Standard Fortran 


Hardware requirements are for the support of tasking and for the provision of adequate 
address space. 
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Features for the Most Popular Tools 

MAC'S SYmbolic MAnipulation (MACSYMA) 

Features: 

1 . Solve equations, factor polynomials, expand functions in series 

2. Equation simplification 

3. Manipulate matrices, trigonometric equation 

4. Integration and differentiation 

5. Automatic code generation for FORTRAN program 

6. File handling system 

7. Interactive help facility 

8. Control function and subroutine procedures 

9. 3D graphics, automatic scaling, and different viewpoints 

10. Interactive sytem 

1 1 . Excellent working environment 

12. Ease of use 

Available computers: 

DEC 10 or 20 
VAX 

Honeywell 6000 Series 
Symbolics LM2 or 3600 

Memory requirements: 

Two Megabytes 
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Symbolic Manipulation Program (SMP) 


Features: 

1 . Solve equations, factor polynomials, expand functions in series 

2. Equation simplification 

3. Manipulate matrices, trigonometric equation 

4. Integration and differentiation 

5. Automatic code generation for FORTRAN or C program 

6. File handling system 

7. Interactive help facility 

8. Control function and subroutine procedures 

9. 3D graphics 

10. Interactive sytem 

1 1 . Excellent working environment 

12. Ease of use 
Available computers: 

VAX 

Memory requirements: 

Two megabytes 

REDUCE II 
Features: 

1 . Solve equations, factor polynomials 

2. Equation simplification 

3. Manipulate matrices 
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4. Integration and differentiation 

5. Automatic code generation for FORTRAN program 

6. File handling system 

7. Control function and subroutine procedure 

8. Interactive system 

Available Computers: 

IBM-360 or 370 
DEC 10 or 20 
VAX 

UNI VAC 1100 
CDC CYBER 
CRAY-1 

BURROUGHS 6700 

APOLLO DOMAIN AND OTHERS 


Memory requirements: 
350 Kilobytes 


MuMath 

Features: 

1. Solve equations, factor polynomials 

2. Equation simplification 

3. Manipulate matrices 

4. Integration and differentiation 

5. Interactive system 
Available computers: 

IBM PC or compatible system 
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VII. MAINTAINING AND VERIFYING 


One of the most difficult problems in software today is the verification and maintenance of 
existing programs, especially programs built up over time with many programmers involved. This 
is particularly true for simulation programs. The traditional method of modifying simulation 
programs by rebuilding and recoding the model is a potential source of many errors. 

The automatic system can eliminate most of these problems. If in error, the only place 
which need be considered to change is the knowledge base at the front, then the rest of the equa- 
tion and program can automatically be changed. The procedure becomes much simpler, which 
makes it easy for the user to maintain and modify the model and program directly. 

Improving maintenance and verification is a major goal for this type of model generation. 
Allowing maintenance to be concerned with equations and logic rather than the implementation is a 
goal for maintenance. 


VIII. GENERAL APPLICATION 


The AI technique for automatic mathematical modeling can apply to many problems which 
are required to be solved: 

1) Complicated mathematical equations derivations that are hard or impossible for humans 
to do manually. 

2) Equations which share the generic theory. 

3) Models which need to be changed considerably. 

Robotics, contact dynamics, and other problems which require a lot of complicated and long 
equation derivation in Matrix Multiplication can very easily be solved symbolically. 

The Space Shuttle Main Engine Simulation contains several equations where there are sets 
of equations sharing the generic theory. The engine model is required to change when the engine 
design in hardware is modified. This is the case used for the Automatic Modeling Example. 

There are many simulation models which need to be changed under different conditions. If 
they are done manually, it will cost a tremendous amount time to modify the model, change the 
program code, and verify the result before it can be used. The AI technique totally improves the 
environment, and the only place which needs to be modified is the knowledge base; the rest of the 
product can be automatically generated. 

Spacelab training applications where several physical models will be developed for a single 
flight. Automatic generation of individual models from equations and logic will eliminate much of 
the repetitive system integration. 
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IX. EXAMPLE APPLICATION 


A demonstration program, called Propulsion System Automatic Modeling (PSAM), was 
designed for modeling the Space Shuttle Main Engine Simulation Mathematic Model. PSAM is 
written in LISP and MACSYMA and runs on a Symbolics 3670 LISP machine. 

The design goals for PSAM were to develop automatic modeling skills for propulsion sys- 
tems and other scientific and engineering applications. The old Engine Model was used as an 
example to study. 

PSAM includes the following features: 

1) User-friendly interface 

2) Automatic Knowledge Base generation 

3) Automatic Equation and Coding generation 

4) Automatic documentation generation. 

The Space Shuttle Main Engine Simulation model was built up from the component process 
elements and their combination into the subprograms. 

The component process elements are pump, hot gas turbine, hydraulic turbine, turbopump, 
combustor, valve, incompressible propellant flow, injector volume with priming for start, hot gas 
heat transfer, and regenerative cooling flow. The subprograms are fuel, oxidizer, and hot gas. 

There are two types of information for a PSAM knowledge base. One is the component 
process elements generic equations, and the other is the information base for the combination of 
the Space Shuttle Main Engine model subprograms and component process elements. 

The system collects the detailed requirment and generates a set of specific equations for the 
component process elements and subprograms. 

PSAM has the ability to: 

1) Create or maintain the knowledge base 

2) Load different knowledge 

3) Automatically generate equations 

4) Output generated equation or FORTRAN code to disk file or option for printout of the 
Laser Printer 

5) Automatic documentation. 
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The FORTRAN code is in generic conventional program format only for the equations. 
Another part of the LISP program will combine all the subprograms and component process 
elements equations and append the header program codes to become a whole compiled program. 
The coding format for the equations are the same for many available conventional software 
languages such as FORTRAN, C, or ADA, etc. So the final program can be built on whichever 
conventional software the user requires by changing only the header requirements and I/O 
functions. 

Several versions and configurations of SSME will be generated and run with a version of 
the engine controller software. 

Another example program was designed for automated kinematic equation generation for 
robot manipulators. Any manipulator can be considered to consist of a series of links connected 
together by joints. Historically, the homogeneous transformation describing the relation between 
one link to the next has been called an A matrix. The position and orientation of the nth link in 
base coordinates are given by the matrix product 


Tn = A1 * A2 * A3 * ... An 


These products of A matrices have been called T matrices. An n-link manipulator can have 
n degrees of freedom, one for each link, and can be positioned and oriented arbitrarily within its 
range of freedom. Most robot controllers equations were developed manually. The complicated 
matrix multiplication is inverse with trigonometric functions. It is tedious and difficult. A small, 
careless mistake can cost major modification of the design. Using the AI technique to derive and 
simplfy the equations symbolically and automatically can reduce a tremendous amount of effort. 

These are the proof of concept examples which lead to future projects. 


X. FUTURE PLANS 


Currently the existing project is being used to test the method. In the future this technique 
can be applied to numerous new projects. The goal is to build the simulation model and maintain it 
automatically. 

The knowledge base will also contain detailed information for automatic documentation. 

The future plans for PSAM are to generate models based upon logic and flow rather than 
equations, and to build control and display modules which will complement the models. 
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APPENDIX 

Sample Problem Solving for 

• Space Shuttle Main Engine Math Model 

• Robot Manipulator 
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SSME Math Model 


Knowledge Base 
Generic Equations 


Information Base 



ORIGINAL PAGE IS 
OF POOR QUALITY 


Generic Equations 


(D105) 

(C106) 


A1 : >PSAM>Macsyma-kb>psam-generic-equat ion. out . 1 
BATCH ( "al : >p3am>macsyma-kb>p9am-generic-equat ion . kb") ; 


(Cl 07) "**********PUMP*************"$ 

(C108) "Pump Flow Variable"? 

(C109) PFV(F,B,DW,S) 

Block ( F - B * (DW /S)); 

(D109) 


B DW 

PFV (F, B, DW, S) J' BLOCK (F = ) 

S 


(Cl 10 ) "Pump Total Pressure"? 

{Clll) PTP (PP ,P,B,S,T) S' 

Block { PP ' P + B * (S) A 2 * T); 

(Dill) 

(Cl 12) "Pump Torque" $ 

(C113) PT (R,B,S,T) :* 

Block ( R - B * ( S ) A 2 * T ) ; 

( D 1 1 3 ) 

(C114 ) "**********TURBINE* *********"? 

(CHS) "Turbine Torque"? 

(C116) TT(R,B,P,T) 

Block { R - B * P * T) ; 

(D116) 

(C117) "Turbine Speed Parameter"? 

(C118) TSP (M,B, S,T) 

Block ( M - B • S / (T A (1 /2) ) ) ? 

(D118) 


2 

PTP (PP, P, B, S, T) := BLOCK (PP = P + B S T) 


2 

PT(R, B, S, T) := BLOCK (R * B S T) 


TT (R, B, P, T) S' BLOCK (R - B P T) 


B S 

TSP (M, B, S, T) BLOCK (M « ) 

1/2 

T 


(C119) "Turbine Weight Flowrate"? 

(C120) TWF (DW,B, dP, r) 

Block ( DW - (B * (dP) • r) A (l/2))? 

(D120 ) TWF ( DW , B, DP , R) 

(C121) "********* * TURBOP UMP SPEED******* "? 


1/2 

BLOCK (DW - (B DP R) ) 


(C122) "Turbo Pump Speed"? 

(C123) TPS <S,B,Tq,S0) s- 

Block ( S - B *' Integrate { Tq,t) + SO); 
(D123) 


/ 

t 

TPS (S, B, TQ, SO) s- BLOCK (S - B I TQ dT + SO) 

] 

/ 


(C124) "**********VALVE** * *******"$ 
(C125) "Valve Area"? 


(C126) VA (A, T, X) 

Block ( A - T * X); 

(D126) 

(Cl 27) -**********FLOW**********"$ 


VA (A, T, X) S' BLOCK (A - T X) 


(C128) "Fuel Flow"? 


(C129) FF (DWf ,Bi, P9, P,B j,R, DWO) 

Block ( DWfsDWf(t), DWO s DWO (t), DWf - B* 'integrate ( { (P9 - P 


(Bi / R)) * DWf A 2),t) + DWO); 

/ 

[ 


BI 
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2 

(D129) FF (DWF, BI, P9, P, BJ, R, DWO) := BLOCK (DWF : DWF(T), DWO : DWO(T), DWF = B I (P9 - P 

) DWF dT + DWO) 

] R 

/ 


(Cl 30 ) "Oxidizer Flow"? 

(C131 ) OF(DWo,Bk,Ppos,P,Bl,A,Ab,Bm, DWi, DWO) := 

Block ( DWo : DWo { t ) , DWi : DWi (t ) , DWO : DWO (t ) , DWo - Bk * 'integrate ( (Ppos - p - Bl * (DWo /(A/Ab)) / '2 - Bm* (DW 
i) *2) , t ) + DWO) ; 

(D1 31 ) OF {DWO, BK, PPOS, P, BL, A, AB, BM, DWI, DWO) := BLOCK (DWO : DWO{T), DWI : DWI (T) , DWO : DWO(T), 


DWO 2 2 

- BL ( ) - BM DWI ) dT + DWO) 

A 

AB 


/ 

t 

DWO = BK I (PPOS - P 
] 

/ 


(C132) "Variable in injector Priming function, Dimensionless" $ 

(Cl 33 ) VIPF (E, Bh, DWO, DWi, EO) : = 

Block (DWO: DWO (t) , DWi:DWi(t), E = Bh * 'integrate ({DWo - DWi),t) + EO) ; 

/ 

[ 

(D133) VIPF (E, BH, DWO, DWI, EO) := BLOCK (DWO : DWO(T), DWI : DWI (T) , E = BH I (DWO - DWI) 

dT + EO) 

1 

/ 


(Cl 34 ) "Weight Flowrate Injector"? 

(C135) WFI (DWi, DWo,Eo) := 

Block { DWi = DWo * Eo) ; 

(D135) WFI (DWI, DWO, EO) := BLOCK(DWI = DWO EO) 

(Cl 3 6) "*★★****★ **COMBUSTER**** ★***★*"$ 

(C137) "Combuster Total Pressure"? 


(C138) CTP (P, Bi , DWI, DW2 , B j , DW3, PO ) : =Block { DWl:DWl(t), DW2:DW2(t), DW3:DW3(t),P = Bi • 'integrate ({ DWI 
+ DW2 - Bj • DW3 ) , t ) + pO); 

/ 

[ 

(D138) CTP (P , BI, DWI, DW2 , BJ, DW3 , PO) : = BLOCK (DWI : DWi (T) , DW2 : DW2 (T) , DW3 : DW3 (T) , P = BI I (DWI 
+ DW2 - BJ DW3 ) dT + PO) 

] 

/ 


(Cl 39) "Combuster Fuel Weight Flowrate"? 

(Cl 40 ) WFF (F, DW2 , DW) : =Block ( F = DW2 / (DWI + DW2) ) ; 

(D140) WFF (F, DW2, DW) 


DW2 

:= BLOCK (F = ) 

DWI + DW2 


(Cl 41 ) "Combuster Temperature"? 

(Cl 42 ) CT (T, TT, Bk, T9) : =Block ( T = TT + Bk * T9) ; 

(D142) CT (T, TT, BK, T9) := BLOCK (T = TT + BK T9) 

(Cl 43 ) "******** ★★COOLING ELEMENT**********"? 


(Cl 44 ) "Cooling Element Total Pressure"? 

(Cl 45) CETP (P, B, DQwl, DQw2, H3, DWi, H, DW, PO) := 

Block { DQwl : DQwl (t ) , DQw2 : DQw2 (t ) , DWi : DWi (t ) , DW: DW { t ) , P = B * ' integrate ( (DQwl + DQw2 + H3 • DWi - H • DW) 
,t) + PO); 

{ D1 45 ) CETP (P , B, DQWl , DQW2 , H3, DWI, H, DW, PO) := BLOCK (DQW1 : DQW1 (T) , DQW2 : DQW2 (T) , DWI : DWI (T) , 
DW : DW (T) , 


/ 

[ 

P = B I (DQW1 


i 


I 
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+ DQW 2 + H3 DWI - H DW) dT + PO) 


ORIGINAL PAGE IS 
OF POOR QUALITY 


i 

/ 


(Cl 46) "Cooling Element Specific Enthalpy"$ 

(C147) CESE(H,B,T) : -Block ( H * B • T) ; 

(D147) CESE (H, B, T) :« BLOCK (H - B T) 

(C148) "Cooling Element Weight Flowrate Main Chamber Heat Exchanger"? 

(Cl 4 9) CEWFM(DW,Bi,P,P5,Bj,R5,DW,DW0) : = 

Block ( DW : DW ( t ) , DW « Bi • 'integrate ( (P - P5 - (Bj / R5) * DW rt 2),t) +DWO); 

/ 2 
[ BJ DW 

(D149) CEWFM (DW, BI, P, P5, BJ, R5, DW, DWO) := BLOCK (DW : DW(T), DW * BI I (P - P5 

) dT + DWO) 

) R5 

/ 

(C150) "Cooling Element Density"? 

(C151 ) CED (R, B, DWi , DW, RO) : = 

Block ( DWisDWi(t), DW:DW(t),R * B * 'integrate ((DWi - DW),t) + RO) ; 

/ 

[ 

(D151) CED (R, B, DWI, DW, RO) :* BLOCK (DWI : DWI (T) , DW : DW(T), R * B I (DWI - DW) dT 

+ RO) 

] 

/ 

(C152) "Cooling Element Temperature"? 

(C153) CET (T, B, P, R) :«* 

Block ( T « B * ( P / R)>; 

B P 

(D153) CET (T, B, P, R) := BLOCK (T ) 

R 

(C154) "Cooling Element Heat Transfer Rate Hot Gas Wall"? 

(C155) HT (DQw, B, T, Tw, DWi ) : = 

Block ( DQw * B * (1.0 + 0.002 * T) * (Tw-T) * (DWi) *0.8); 

0.8 

(D155) HT (DQW, B, T, TW, DWI) :« BLOCK (DQW * B (1.0 + 0.002 T) (TW - T) DWI ) 

(C156) "Cooling Element Heat Transfer Rate TC"? 

(C157 ) HTTC (DQtc, B, Tc, Twl, DWcn) 

Block ( DQtc - B * (Tc - Twl) * DWcn*0.8); 

0.8 

(D157 ) HTTC (DQTC, B, TC, TW1 , DWCN) BLOCK (DQTC * B (TC - TW1) DWCN ) 

(C158) "Cooling Element Hot Gas Wall Temperature"? 

(Cl 59 ) HGT (TT,B, DQtc, DQwl, TO) : - 

Block ( DQtc : DQtc (t) , DQwl : DQwl (t ) , TT ** B • 'integrate ((DQtc - DQwl) ,t) + TO); 

/ 

[ 

(D159) HGT (TT, B, DQTC, DQW1, TO) BLOCK (DQTC : DQTC (T) , DQW1 : DQW1 (T) , TT *= B I (DQTC - D 

QW1 ) dT + TO) 

] 

/ 

(Cl 60) "Cooling Element Ambient thrust Chamber Temperature"? 

(C161 ) ATT (TT, B, DQw2 , TO) : “Block ( DQw2 : DQw2 (t ) , TT =* B * 'integrate ((-DQw2),t> +T0) ; 

/ 

c 

(D161) ATT (TT, B, DQW2, TO) : = BLOCK (DQW2 : DQW2 (T) , TT = 8 I (- DQW2) dT + TO) 

1 

/ 

(Cl 62) ****** *****Flowmeter*** ********$ 

(C163) "Volumetric Flowrate"? 

(Cl 64) VF(Q,C1,T,C2,C3,C4,DW,C5,S) : = 

Block ( Q * [1.0*- Cl *( T - C2 )]*( C3 + C4 • DW + C5 * SJ); 

(D164) VF(Q, Cl, T, C2, C3, C4, DW, C5, S) BLOCK (Q = [1.0 - Cl (T - C2)J [C3 + C4 DW 
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+ C5 S] ) 

(Cl 65 ) "Total Pressure"? 

(Cl 66 ) TP (P,Pf,RF,DW,R) : = 

Block ( P = Pf - (RF * DVT2)/R); 

2 

RF DW 

(D166) TP (P, PF, RF, DW, R) := BLOCK (P - PF - ) 

R 

(Cl 67 ) "Flowmeter Temperature"? 

(Cl 68 ) FT (T, Ti , S, Coe ) : = 

Block ( T = Ti + Coe * S) ; 

( D1 68 ) FT (T, TI, S, COE) := BLOCK (T = TI + COE S) 

(Cl 69 ) ********** * ***p re burner Fuel Supply Duct*****************"? 

(C170) "Supply Duct Total Pressure"? 

(C171 ) SDTP (P, B, DW1 , DW2 , DW3 , DW4 , PO) := 

Block ( dwl:dwl(t), dw2:dw2(t), dw3:dw3(t), dw4:dw4(t),P = B • ' integrate < (DWl - DW2 - DW3 +DW4),t) + PO); 
(D171) SDTP (P, B, DWl, DW2 , DW3, DW4 , PO) := BLOCK (DWl : DWl (T) , DW2 : DW2(T), DW3 : DW3 (T) , DW4 : DW4 (T) 


(DWl - DW2 - DW3 + DW4 ) dT + PO) 

(C172) "Supply Duct Temperature"? 

(Cl 73) SDT (T, B, P, R) : = 

Block ( T = B * (P / R) ) ; 

(D173) 


/ 

[ 

P = B I 
] 

/ 


B P 

SDT (T, B, P, R) := BLOCK (T = ) 

R 


(C174 ) "Supply Duct Density"? 

(C175) SDD (R, DWl, DW2, Rr, B):= 

Block ( R - (DWl + DW2) / ( (DW2/Rr) + B * DWl]); 

[DWl + DW2 ) 

(D175 ) SDD (R, DWl, DW2, RR, B) := BLOCK (R ) 

DW2 

[ + B DWl] 

RR 


(C176) "****************Hi gh Pressure Oxidizer Pump Preburner Boost stage**************"? 

(C177) "High Pressure Oxidizer boostpump Flow variable"? 

(Cl 78) HPOPFV(F,B,DW,S) : = 

Block (F =» B * (DW / S) ) ; 

B DW 

(D178) HPOPFV (F, B, DW, S) := BLOCK (F = ) 

S 


(C179) "High Pressure Oxidizer Boostpump Torque"? 

(C180) HPOPT (T, B, S, TT, F) : = 

Block ( T - B44 * (S^2) * TT * (F) ) ; 

2 

(D180) HPOPT (T, B, S, TT, F) := BLOCK (T = B44 S TT F) 


(Cl 81 ) "High Pressure Oxidizer Boostpump Pressure"? 

(Cl 82 ) HPOPP (P3 , B, S, TT, F, P2 ) : * 

Block ( P3 - B * ( S*2) • TT • F + P2); 

2 

(D182) HPOPP (P3, B, S, TT, F, P2) := BLOCK (P3 = B S TT F + P2) 

(C183) "High Pressure Oxidizer Boostpump Weight Flowrate"? 

(C184) HPOPW (DW3, Bl, P3, P, B2, DWO) 

Block ( DW3 : DW3 (t ) , DW3 = Bl • 'integrate ( (P3 - P - B2* DW3*DW3) , t) + DWO ); 


/ 
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t 

(D1 04 ) HPOPW (DW3, Bl, P3, P, B2 , DWO) := BLOCK(DW3 : DW3 (T) , DW3 e B1 I {P3 - P - B2 DW3 DW3 

) dT + DWO> 

] 

/ 


(C185 ) ******************* *Main Chamber Injector, Fuel and Hot Gas *************************«$ 

(Cl 86) "Main Chamber Injector Total Pressure"? 

(C187) MCITP (P, B1 , DHl , DW2, DW3, B2, DW4, P0): = 

Block (DW1:DW1 (t) , DW2:DW2(t), DW3:DW3(t), DW4:DW4(t), P = B1 • ' integrate ( (DW1 + DW2 + DW3 - B2* DW4),t) 
+ PO) ; 

(D187) MCITP (P, Bl, DW1, DW2, DW3, B2, DW4, PO) := BLOCK (DW1 : DW1 <T) , DW2 : DW2 (T) , DW3 : DW3 (T) , DW4 : 
DW4 (T) , 


1 + DW2 + DW3 - B2 DW4 ) dT + PO) 


/ 

[ 

P = Bl I (DW 

J 

/ 


(Cl 88 ) "Main Chamber Injector Weight Flowrate"? 

(Cl 89 ) MCIW(DW,B,P1,T,PC,P2,TT) :« 

Block (DW - B * PI * ( T * PC/P2) /TT^O. 5) ; 

B PI (T PC) 


(D189) 


MCIW (DW, B, PI, T, PC, P2, TT) := BLOCK (DW 


P2 

) 

0.5 

TT 


(C190) "Main Chamber Injector Temperature"? 

(Cl 91 ) MCIT (TFI, Bl , TFP, B2, TOP, B3, T5) : = 

Block (TFI - Bl * TFP + B2 • TOP + B3 • T5) ; 

(D191) MCIT (TFI, Bl, TFP, B2, TOP, B3, T5) :« BLOCK (TFI - Bl TFP + B2 TOP + B3 T5 

) 

(C192) "********************Main Chamber Injector, Oxidizer*********************************"? 

(Cl 93) "Main Chamber Injector, Oxidizer Weight Flowrate"? 

(Cl 94) MCWF (DWoi, DW, E) : * 

Block (DWoi - DW • E); 

(D194) MCWF (DWOI, DW, E> := BLOCK (DWOI - DW E) 

(C195) "Main Chamber Ir jector. Oxidizer Weight"? 

(Cl 96) MCW (W, DW1 , DW2 , WO):- 

Block (DW1:DW1 (t) , DW2:DW2(t), W - ' integrate ( (DW1 - DW2),t) + WO) ; 

/ 

I 

(D196) MCW (W, DW1 , DW2, WO) : * BLOCK (DW1 : DW1 (T) , DW2 : DW2(T), W = I (DW1 - DW2) dT 

+ WO) 

1 

/ 


(Cl 97) "Main Chamber injector. Oxidizer Variable in injector Priming function"? 

(Cl 98) MCV (E, T, W) : = 

Block (E - T * W); 

(D198) MCV ( E , T, W) := BLOCK (E - T W) 

(C199) «************************Min Thrust chamber*******************************************"? 

(C200) "Combuster Injector end static Total Pressure"? 

(C201) CITP (P, B, P) : * 

Block (P - B * P) f 

(D201) CITP (P, B, P) :« BLOCK (P -BP) 

(C202) "Main thrust Chamber Total Pressure"? 

(C203) MTCTP (P, Bl, B2, DW1, DW2, DW3, PO) :« 

Block (DW1:DW1 (t) , DW2:DW2(t), DW3:DW3(t), P - Bl * * integrate ( (B2*DW1 + DW2 - DW3),t) + PO) ; 
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(D203) MTCTP {P, Bl, B2, DW1, DW2, DW3, PO) := BLOCK (DW1 : DW1 (T) , DW2 : DW2 (T) , DW3 ; DW3 (T ) , P = Bl I (B 
2 DW1 + DW2 - DW3) dT + PO) 

1 

/ 


(C204) "A******************* ***Nozzle Bypass Element***************************************"? 
(C205) "Nozzle Heat exchanger bypass Weight Flowrate"? 

(C206) NHWF (DW, Bl , P, P9, B2, DW1, A, DWO) := 

Block (DW1 : DW1 (t ) , DW = B1 * # integrate { (P - P9 - B2 * (DW1/ (A) A 2) ) , t) + DWO); 

/ 

[ 

1 

(D206) NHWF (DW, Bl , P, P9, B2, DW1, A, DWO) := BLOCK (DW1 : DW1 (T) , DW = B1 I <P - P9 

-) dT + DWO) 

] 

/ 


B2 DW 


2 


A 


(C207) "Coolant Control Valve"? 

(C208) CCV (DW , DW1) : = 

Block {DW = DW1 ) ; 

{D208 ) CCV (DW, DW1 ) := BT,OCK(DW = DW1 ) 

(D209) DO»\t 

(C210) CLOSEFILEO; 
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I 

(Deframe Fop3 

(Class BOOST-STAGE) 

(Unit "Flow-Variable") 

(Nomenclature "High Pressure Oxidizer boostpump Flow variable") 

(Input -parameter B43 DWop3 So2) ) 
j (Deframe Top3 

, (Class BOOST-STAGE) 

(Unit "Torque") 

(Nomenclature "High Pressure Oxidizer Boostpump Torque") 

(Input-parameter B44 So2 Trop3 Fop3) ) 

(Deframe Pod3 

(Class BOOST-STAGE) 

(Unit "Pressure") 

(Nomenclature "High Pressure Oxidizer Boostpump Pressure") 

(Input -parameter B45 So2 Tpop3 Fop3 Pod2) ) 

(Deframe DWop3 

(Class BOOST-STAGE) 

(Unit "Weight-Flowrate") 

(Nomenclature "High Pressure Oxidizer Boostpump Weight Flowrate") 
(Input-parameter B46 Pod3 Ppos B104 DWop30)) 

(Deframe Trop3 

(Class BOOST-STAGE) 

(Unit "Table-Function") 

(Nomenclature "High Pressure Oxidizer Boostpump Torque") 

(Input -parameter )) 

(Deframe Tpop3 

(Class BOOST-STAGE) 

(Unit "Table-Function") 

(Nomenclature "high Pressure Oxidizer boostpump Pressure Rise Character") 
(Input-parameter )) 

(Deframe DWop30 

(Class BOOST-STAGE) 

(Unit "Initial-Condition") 

(Nomenclature "DWop3 when t=0") 

(Input-parameter 105.43)) 
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{Deframe Pfp 

(Class COMBUSTER) 

(Unit "Combust er- Total -Pres sure" ) 

(Nomenclature "Fuel Preburner Total Pressure") 

(Input-parameter B56 DWfpf DWfpoi Bill DWft2 PfpO) ) 

(Deframe Pop 

(Class COMBUSTER) 

(Unit " Combu st er-Total -Pres sure" ) 

(Nomenclature "Oxidizer Preburner Total Pressure") 
(Input-parameter b63 DWopf DWopoi B107 DWot2 PopO) ) 

(Deframe Ffp 

(Class COMBUSTER) 

(Unit "Weight -Flowrate-Fraction" ) 

(Nomenclature "Fuel Preburner Flowrate Fraction") 
(Input-parameter DWfpoi DWfpoi+DWfpf ) ) 

(Deframe Fop 

(Class COMBUSTER) 

(Unit "Weight-Flowrate-Fract ion" ) 

(Nomenclature "Oxidizer Preburner Flowrate Fraction") 
(Input-parameter DWopoi DWopoi+DWopf ) ) 

(Deframe Tfp 

(Class COMBUSTER) 

(Unit "Combuster-Temperature") 

(Nomenclature "Fuel Preburner Combuster Temperature”) 
(Input-parameter TTpfp B112 T9) ) 

(Deframe Top 

(Class COMBUSTER) 

(Unit "Combuster-Temperature") 

(Nomenclature "Oxidizer Preburner Combuster Temperature") 

(Input -parameter TTfop B113 T9) ) 

(Deframe PfpO 

(Class COMBUSTER) 

(Unit "Initial-Condition") 

(Nomenclature "Pfp when t=0") 

(Input-parameter 5940.0)) 

(Deframe PopO 

(Class COMBUSTER) 

(Unit "Initial-Condition") 

(Nomenclature "Pop when t=0") 

(Input-parameter 5929.9)) 

(Deframe TTffp 

(Class COMBUSTER) 

(Unit "Table-Function") 

(Nomenclature "Fuel Preburner Hot Gass Temperature Function") 
(Input-parameter ) ) 

(Deframe TTfop 

(Class COMBUSTER) 

(Unit "Table-Function”) 

(Nomenclature "Oxidizer Preburner Hot Gas Temperature Function") 
(Input-parameter ) ) 
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{Deframe 


(Deframe 


(Deframe 


(Deframe 


(Deframe 


(Deframe 


(Deframe 


(Deframe 


(Deframe 


(De frame 


(Deframe 


(Deframe 


(Deframe 


(Deframe 


(Deframe 


(De frame 
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P5 

(Class COOLING) 

(Unit "Total-Pressure") 

(Nomenclature "Main Chamber Heat Exchanger Total Pressure") 
(Input-parameter 1 DQwl5 DQw25 H3 DWmc H5 DW5 P50) ) 

PA 

(Class COOLING) 

(Unit "Total-Pressure") 

(Nomenclature "Nozzle Regenerative Exchanger Element Total Pressure") 
(Input -parameter b84 DQwl4 DQw24 H4 DW4 H3 DWfn P40) ) 


H5 

(Class COOLING) 

(Unit "Specific-Enthalpy") 

(Nomenclature "Main Chamber Heat Exchanger Specific Enthalpy") 
(Input-parameter blOl T5) ) 

H4 

(Class COOLING) 

(Unit "Specific-Enthalpy") 

(Nomenclature "Nozzle Regenerative Cooling Element Specific Enthalpy") 
(Input-parameter B102 T4) ) 


DW5 

(Class COOLING) 

(Unit "Weight-Flowrate") 

(Nomenclature "Main Chamber Heat Weight Flowrate") 
(Input-parameter B75 Pmfvd P5 B76 R5 DWmc DWmcO) ) 

DW4 

(Class COOLING) 

(Unit "Weight -Flowrate") 

(Nomenclature "Nozzle Regenerative Cooling Weight Flowrate") 
(Input -parameter B94 P4 P9 B95 R4 DW4 DW40)) 

R5 

(Class COOLING) 

(Unit "Density") 

(Nomenclature "Main Chamber Element Density") 
(Input-parameter B77 DWmc DW5 R50) ) 


R4 

(Class COOLING) 

(Unit "Density") 

(Nomenclature "Nozzle Regenerative Element Density") 

(Input-parameter B87 DWfn DW4 R40) ) 

T5 

(Class COOLING) 

(Unit "Temperature") 

(Nomenclature "Main Chamber Heat Exchanger Temperature") 
(Input-parameter B78 P5 R5) ) 

T4 

(Class COOLING) 

(Unit "Temperature") 

(Nomenclature "Nozzle Regenerative Cooling Temperature") 

{ Input -parameter B88 P4 R4) ) 

DQwl5 

(Class COOLING) 

(Unit "Heat-Transfer") 

(Nomenclature "Main Chamber Hot Gas Heat Transfer Rate") 
(Input-parameter B79 T5 TW15 DWmc) ) 

DQw25 

(Class COOLING) 

(Unit "Heat-Transfer") 

(Nomenclature "Main Chamber Ambient Thrust Chamber Heat Transfer Rate") 
(Input -parameter B80 T5 TW25 DWmc) ) 

DQwl4 

(Class COOLING) 

(Unit "Heat-Transfer") 

(Nomenclature "Nozzle Regenerative Hot Gas Heat Transfer Rate") 
(Input-parameter B89 T4 TW14 DWfn) ) 

DQw24 

(Class COOLING) 

(Unit "Heat-Transfer") 

(Nomenclature "Nozzle Regenrative Ambient Thrust Heat Transfer Rate") 
(Input-parameter B90 T4 TW24DWfn)) 

DQtc5 

(Class COOLING) 

(Unit "Heat-Transfer-TC") 

(Nomenclature "Main Chamber Cooling Element Heat Transfer Rate TC") 
(Input-parameter B81 Tc Twl5 DWcn) ) 

DQtc4 

(Class COOLING) 

(Unit "Heat-Transfer-TC") 

(Nomenclature "Nozzle Regenerative Heat Transfer Rate TC") 



(Input -parameter B91 Tc Twl4 DWcn) ) 

(Deframe Twl5 

(Class COOLING) 

(Unit "Hot-Gas-Temperature" ) 

(Nomenclature "Hot Gas Temperature") 
(Input-parameter B82 DQtcS DQwl5 Twl50) ) 

(Deframe Tw25 

(Class COOLING) 

(Unit "Ambient-Thrust-Temperature” ) 

(Nomenclature "Ambient thrust chamber Temperature") 
(Input-parameter B83 DQw25 Tw250) ) 

(De frame Twl4 

(Class COOLING) 

(Unit "Hot-Gas-Temperature") 

(Nomenclature "Hot Gas temperature") 
(Input-parameter B92 DQtc4 DQwl4 Twl40) ) 

(Deframe Tw24 

(Class COOLING) 

(Unit "Ambient -Thrust -Temperature") 

(Nomenclature "Ambient Thrust Temperature ”) 
(Input-parameter B93 DQw24 Tw240)) 

(Deframe P50 

(Class COOLING) 

(Unit "Initial-Condition") 

(Nomenclature "P5 when t=0") 

(Input-parameter 5143.0)) 

(Deframe P40 

(Class COOLING) 

(Unit "Initial-Condition") 

(Nomenclature "P4 when t=0") 

(Input-parameter 6794.0)) 

(De frame DWmcO 

(Class COOLING) 

(Unit "Initial-Condition") 

(Nomenclature "DWmc when t=0") 

(Input-parameter 35.24)) 

(De frame DW40 

(Class COOLING) 

(Unit "Initial-Condition") 

(Nomenclature "DW4 when t=0") 

(Input-parameter 41.98)) 

(Deframe R50 

(Class COOLING) 

(Unit "Initial-Condition") 

(Nomenclature "R5 when t=0") 

(Input-parameter 0.001402)) 

(Deframe R40 

(Class COOLING) 

(Unit "Initial-Condition") 

(Nomenclature "R4 when t=0") 

(Input-parameter 6650.0)) 

(Deframe Twl50 

(Class COOLING) 

(Unit "Initial-Condition") 

(Nomenclature "TW15 when t=0") 

(Input-parameter 1330.0)) 

(Deframe Tw250 

(Class COOLING) 

(Unit "Initial-Condition") 

(Nomenclature "Tw25 when t=0") 

(Input-parameter 544.5)) 

(Deframe Twl40 

(Class COOLING) 

(Unit "Initial-Condition”) 

(Nomenclature "Twl4 when t=0") 

(Input-parameter 1260.0)) 

(Deframe Tw240 

(Class COOLING) 

(Unit "Initial-Condition") 

(Nomenclature "Tw24 when t=0") 

(Input-parameter 679.3)) 


ORIGINAL PAGE IS 
OF POOR QUALITY 


24 


(Deframe 


(De frame 


(De frame 


(Deframe 


(De frame 


(Deframe 


(De frame 


) ) 

(Deframe 


(Deframe 


(Deframe 


(Deframe 


(Deframe 


(Deframe 


(Deframe 
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DWfpf 

(Class FLOW) 

(Unit "Fuel-Flow") 

(Nomenclature "Fuel Preburner Fuel Flow") 

(Input -parameter B50 P9 Pfp B51 R9 DWfpfO) ) 

DWopf 

(Class FLOW) 

(Unit "Fuel -Flow") 

(Nomenclature "Oxidizer Preburner Fuel Flow") 

(Input -parameter B57 P9 Pop B58 R9 DWopfO)) 

DWfpo 

(Class FLOW) 

(Unit "Oxidizer-Flow") 

(Nomenclature "Fuel Preburner Oxidizer Weight Flowrate") 

(Input -parameter B52 Ppos Pfp B53 A Ab B54 DWfpoi DWfpoO) ) 

DWopo 

(Class FLOW) 

(Unit "Oxidizer-Flow") 

(Nomenclature "Oxidizer Preburner Oxidizer Weight Flowrate") 
(Input-Parameter B59 Ppos Pop B61 A Ab B60 DWopoi DWopoO) ) 

Efpo 

(Class FLOW) 

(Unit "Variable- in-injector-Priming") 

(Nomenclature "Fuel Preburner Oxidizer Variable in injhector Priming") 
(Input-parameter B55 DWfpo DWfpoi EfpoO) ) 

Eopo 

(Class FLOW) 

(Unit "Variable-in-in jector-P riming" ) 

(Nomenclature "Oxidizer Preburner Oxidizer Variable in injector Priming") 
(Input -parameter B62 DWopo DWopoi EopoO) ) 

DWfpoi 
(Class FLOW) 

(Unit "Weight -Flowrate-In jector") 

(Nomenclature "Fuel Preburner Oxidizer Injector Weight Flowrate") 
(Input-parameter DWfpo Efpo 


DWopoi 
(Class FLOW) 

(Unit "Weight-Flowrate -In jector" ) 

(Nomenclature "Oxidizer Preburner Oxidizer Injector Weight Flowrate") 
(Input-parameter DWopo Eopo) ) 

DWfpfO 
(Class PUMP) 

(Unit "Initial-Condition") 

(Nomenclature "DWfpf when t-0") 

(Input -parameter 83.15)) 

DWopfO 
(Class PUMP) 

(Unit "Initial-Condition") 

(Nomenclature "DWopf when t=0") 

(Input -parameter 38.93)) 

DWfpoO 

(Class PUMP) 

(Unit "Initial-Condition") 

(Nomenclature "DWfpo when t=0") 

(Input-parameter 79.30)) 

DWopoO 
(Class PUMP) 

(Unit "Initial-Condition") 

(Nomenclature "DWopo when t=0") 

(Input-parameter 26.13)) 

EfpoO 

(Class PUMP) 

(Unit "Initial-Condition" ) 

(Nomenclature "EFpo when t=0") 

(Input-parameter 1.0)) 

EopoO 

(Class PUMP) 

(Unit "Initial-Condition") 

(Nomenclature "Eopo when t=0") 

(Input-parameter 1.0)) 



(Deframe QF 

(Class FLOWMETER) 

(Unit "Volumetric-Flowrate") 

(Nomenclature "Fuel Volumetric Flowrate") 

(Input-parameter 0.008783 Tfs 37.0 -204.852 101.323 DWfd2 0.20206 Sfl) ) 

(Deframe Pffm 

(Class FLOWMETER) 

(Unit "Total-Pressure") 

(Nomenclature "Fuel Flowmeter Total Pressure") 

(Input-parameter Pfdl RFfcd DWfd2 Rfdl) ) 

(Deframe Tfpl 

(Class FLOWMETER) 

(Unit "Flowmeter-Temperature") 

(Nomenclature "Low Pressure Fuel Pump Flowmeter Temperature") 

(Input -parameter Tfs 0.0015789 Sfl)) 

(Deframe Tod2 

(Class FLOWMETER) 

(Unit "Flowmeter-Temperature") 

(Nomenclature "Low pressure Oxidizer pump Flowmeter Temperature") 
(Input-parameter Tos 0.00547 So2)) 
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ORIGINAL PA ££ S3 
Op POOR QUALITY 

(Deframe Pfi 

(Class MAIN -CHAMBER) 

(Unit "Pressure" ) 

(Nomenclature "Main Chamber Injector, fuel and Hot Gas Total Pressure") 
(Input-parameter B64 DWftl DWot2 DWft2 B108 DWfi PfiO) ) 

(Deframe DWfi 

(Class MAIN-CHAMBER) 

(Unit "Weight-Flowrate”) 

(Nomenclature "Main Chamber injector weight flowrate") 

(Input-parameter B64 Pfi Tpr Pc Ppi Tfi)) 

{Deframe Tfi 

(Class MAIN-CHAMBER) 

(Unit "Temperature" ) 

(Nomenclature "Main Chamber Temperature") 

(Input-parameter B66 Tfp B67 Top B68 T5 ) ) 

(Deframe DWoi 

(Class MAIN-CHAMBER) 

(Unit "Oxidizer-Weight -Flowrate" ) 

(Nomenclature "Main Chamber injector. Oxidizer Weight Flowrate") 

(Input -parameter DWmov Eoi)) 

(Deframe Woi 

(Class MAIN-CHAMBER) 

(Unit "Oxidizer-Weight") 

(Nomenclature "Main Chamber Injector, Oxidizer Weight") 

(Input -parameter DWmov DWoi Woi)) 

(Deframe Eoi 

(Class MAIN-CHAMBER) 

(Unit "Oxidizer-Variable") 

(Nomenclature "Main Chamber injector. Oxidizer variable in injector priming function") 
(Input-parameter TEoi Woi) ) 

(Deframe PfiO 

(Class MAIN-CHAMBER) 

(Unit "Initial -Condition") 

(Nomenclature "Pfi when t=0") 

(Input-parameter 3662.0)) 



(Deframe Pcies 

(Class MAIN-THRUST-CHAMBER) 

(Unit "Combust er-Total -Pres sure” ) 

(Nomenclature "Combuster Injector end static Total Pressure") 
(Input-parameter B114 Pc)) 

(Deframe Pc 

(Class MAIN-THRUST-CHAMBER) 

(Unit "Main-thrust-Total -Pres sure" ) 

(Nomenclature "Main thrust chamber Total Pressure") 
(Input-parameter B72 B109 DWfi DWoi DWcn PcO) ) 

(Deframe PcO 

(Class MAIN-THRUST-CHAMBER) 

(Unit "Initial-Condition") 

(Nomenclature "Pc when t=0") 

(Input -parameter 3330.0)) 
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(Deframe DWfnbp 

(Class NOZZLE-BYPASS) 

(Unit "Nozzle-Heat-Weight-Flowrate") 

(Nomenclature "Nozzle Heat Exchanger bypass weight flowrate") 
(Input -parameter B98 Pmfvd P9 B99 DWfnbp Accv DWfnbpO) ) 
(Deframe DWccv 

(Class NOZZLE-BYPASS) 

(Unit "Coolant-Control-Weight -Flowrate") 

(Nomenclature "Coolant Control Valve”) 

(Input -parameter DWfnbp) ) 

(Deframe DWfnbpO 

(Class NOZZLE-BYPASS) 

(Unit "Initial -Condition") 

(Nomenclature "DWfnbp when t=0") 

(Input -parameter 80.19) ) 
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(Deframe P9 

(Class PREBURNER-SUPPLY-DUCT) 

(Unit "Supply-Duct-Total -Pres sure" ) 

(Nomenclature "Preburner Fuel Supply Element Total Pressure") 
(Input -parameter B96 DW4 DWopf DWfpf DWfnbp P90)) 

(Deframe Ppos 

(Class PREBURNER-SUPPLY-DUCT) 

(Unit "Supply-Duct-Total -Pres sure") 

(Nomenclature "Preburner Oxidizer Supply Total Pressure") 
(Input-parameter B49 DWop3 DWfpo Dwopo DW PposO) ) 

(Deframe T9 

(Class PREBURNER-SUPPLY-DUCT) 

(Unit "Supply-Duct-Temperature" ) 

(Nomenclature "Preburner Fuel Supply Element Temperature") 
(Input-parameter B100 P9 R9) ) 

(Deframe R9 

(Class PREBURNER-SUPPLY-DUCT) 

(Unit "Supply-Duct-Density") 

(Nomenclature "Preburner Fuel Supply Element Density") 
(Input-parameter DWfnbp DW4 R4 BIOS)) 

(Deframe P90 

(Class PREBURNER-SUPPLY-DUCT) 

(Unit "Initial-Condition") 

(Nomenclature "P9 when t=0") 

(Input-parameter 6650.0) ) 

(Deframe PposO 

(Class PREBURNER-SUPPLY-DUCT) 

(Unit "Initial-Condition") 

(Nomenclature "Ppos when t=0”) 

(Input -parameter 8457.0)) 
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(Deframe 


(Deframe 


(Deframe 


(Deframe 


(Deframe 


(Deframe 


(Deframe 


(Deframe 


(Deframe 


(Deframe 


(Deframe 


CPJ C’N/'L P/>Q£ !$ 

OF POOR QUALITY 

Ffpl 

(Class PUMP) 

(Unit "Flow-Variable") 

(Nomenclature "Low Pressure Fuel Pump Flow variable”) 
(Input-parameter Bll DWfdl Sfl) ) 

Ffp2 

(Class PUMP) 

(Unit "Flow-Variable”) 

(Nomenclature "High Pressure Fuel Pump Flow Variable") 
(Input-parameter B18 DWfd2 Sf2) ) 

Fdpl 

(Class PUMP) 

(Unit "Flow-Variable") 

(Nomenclature "Low Pressure Oxidizer Pump Flow variable") 
(Input-parameter B27 DWmov+DWop3 Sol) ) 

F op2 

(Class PUMP) 

(Unit "Flow-Variable") 

(Nomenclature "High Pressure Oxidizer Pump Flow variable") 
(Input-parameter B39 DWmov+DWotl+DWop3 So2) ) 

Pfdl 

(Class PUMP) 

(Unit "Total-Pressure") 

(Nomenclature "Low Pressure Fuel Pump Discharge Total pressure") 
(Input-parameter Pfs B12 Sfl Tpfpl) ) 

Pfd2 

(Class PUMP) 

(Unit "Total-Pressure") 

(Nomenclature "High Pressure Fuel Pump Total Pressure") 
(Input-parameter Pfdl B19 Sf2 Tpfp2) ) 

Podl 

(Class PUMP) 

(Unit "Total-Pressure") 

(Nomenclature "Low Pressure Oxidizer Pump Total Pressure") 
(Input-parameter Pos B28 SOI Tpopl)) 

Pod2 

(Class PUMP) 

(Unit "Total-Pressure") 

(Nomenclature "High Pressure Oxidizer Pump Total Pressure”) 
(Input -parameter Podl B40 So2 Tpop2) ) 

Rfpl 

(Class PUMP) 

(Unit "Torque") 

(Nomenclature "Low Pressure Fuel Pump "Torque"") 

(Input-parameter B13 Sfl Trfpl)) 

Rfp2 

(Class PUMP) 

(Unit "Torque") 

(Nomenclature "High Pressure Fuel Pump "Torque”") 
(Input-parameter B20 Sf2 Trfp2) ) 

Ropl 

(Class PUMP) 

(Unit "Torque") 

(Nomenclature "Low Pressure Oxidizer Pump "Torque"") 

(Input -parameter B34 Sol Tropl)) 


(Deframe Tpfpl * 

(Class PUMP) 

(Unit "Table-Function") 

(Nomenclature "Low Pressure Fuel Pump Pump Pressure Rise Characteristics") 
(Input -parameter )) 

(Deframe Tpfp2 

(Class PUMP) 

(Unit "Table-Function") 

(Nomenclature "High Pressure Fuel Pump Pump Pressure Rise Characteristics") 
(Input-parameter )) 

(Deframe Tpopl 

(Class PUMP) 

(Unit "Table-Function") 

(Nomenclature "Low Pressure Oxidizer Pump") 

(Input-parameter )) 

(Deframe Tpop2 

(Class PUMP) 

(Unit "Table-Function") 

(Nomenclature "High Pressure Oxidizer Pump Pump Pressure Rise Characteristics") 
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(Input -parameter ) ) 

(Deframe Trfpl 

(Class PUMP) 

(Unit "Table-Function”) 

(Nomenclature "Low Pressure Fuel Pump Pump Torque Characteristics") 
(Input-parameter )) 

(Deframe Trfp2 

(Class PUMP) 

(Unit "Table-Function") 

(Nomenclature "High Pressure Fuel Pump Pump Pump Torque Characteristics") 
(Input -parameter ) ) 

(Deframe Tropl 

(Class PUMP) 

(Unit "Table-Function") 

(Nomenclature "Low Pressure Oxidizer Pump Pump Torque Characteristics") 

{ Input -parameter )) 

(Deframe Torque 

(Class PUMP) 

(Unit "Torque") 

(Nomenclature "High Pressure Oxidizer Pump \"Torque\"" ) 

(Input -parameter B41 So2 Trop2)) 

(Deframe Trop2 

(Class PUMP) 

(Unit "Table-Function") 

(Nomenclature "High Pressure Oxidizer Pump Pump Torque Characteristics") 
(Input -parameter )) 


ORIGINAL PAGE IS 
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ORIGINAL PASS 

{Deframe TQftl OF POOR QUALITY 

( Class TURBINE) 

(Unit "Torque”) 

(Nomenclature "Low Pressure Fuel Turbine Torque") 

(Input-parameter B16 P5 Trftl) ) 

(Deframe TQft2 

(Class TURBINE) 

(Unit "Torque") 

(Nomenclature "High Pressure Fuel Turbine Torque") 

(Input-parameter B26 Pfp Trft2) ) 

(Deframe TQotl 

(Class TURBINE) 

(Unit "Torque") 

(Nomenclature "Low Pressure Oxidizer Turbine Torque") 

(Input-parameter B36 DWotl^2 Trotl)) 

(Deframe TQot2 

(Class TURBINE) 

(Unit "Torque") 

(Nomenclature "High Pressure Oxidizer Turbine") 

(Input-parameter B29 Pop Trot2) ) 

(Deframe Uftl 

(Class TURBINE) 

(Unit "Speed-Parameter") 

(Nomenclature "Low Pressure Fuel Turbine Speed Parameter") 

(Input-parameter B15 Sfl T5) ) 

(Deframe Uft2 

(Class TURBINE) 

(Unit "Speed-Parameter") 

(Nomenclature "High Pressure Fuel Turbine Speed Parameter") 

(Input -parameter B25 Sf2 Tfp) ) 

(Deframe Fotl 

(Class TURBINE) 

(Unit "Speed-Parameter") 

(Nomenclature "Low Pressure Oxidizer Turbine Flow Variable") 

(Input-parameter B38 Sol DWotl)) 

(Deframe Uot2 

(Class TURBINE) 

(Unit "Speed-Parameter") 

(Nomenclature "High Pressure Oxidizer Turbine Speed Parameter") 
(Input-parameter B47 So2 Rop*Top) ) 

(Deframe DWftl 

(Class TURBINE) 

(Unit "Weight-Flowrate") 

(Nomenclature "Low Pressure Fuel Turbine Weight Flowrate") 

(Input-parameter B17 P5-Pfi R5) ) 

(Deframe DWft2 

(Class TURBINE) 

(Unit "Weight-Flowrate") 

(Nomenclature "High Pressure Fuel Turbine Weight Flowrate") 

(Input -parameter B24 Pfp Tpr) ) 

(Deframe DWotl 

(Class TURBINE) 

(Unit "Weight-Flowrate") 

(Nomenclature "Low Pressure Oxidizer Turbine Weight Flowrate") 

( Input -parameter <1/(B37 + Rotl)) Pod2-Podl 1)) 

(Deframe DWot2 

(Class TURBINE) 

(Unit "Weight-Flowrate" ) 

(Nomenclature "High Pressure Oxidizer Turbine Weight Flowrate") 
(Input-parameter B48 Pop Tpr/Top) ) 

(De frame Trftl 

(Class TURBINE) 

(Unit "Table-Function") 

(Nomenclature "Low Pressure Fuel Turbine Turbine Torque Characteristics") 
(Input-parameter ) ) 

(Deframe Trft2 

(Class TURBINE) 

(Unit "Table-Function" ) 

(Nomenclature "High Pressure Fuel Turbine Turbine Torque Characteristics") 
(Input-parameter ) ) 

(Deframe Trotl 

(Class TURBINE) 

(Unit "Table-Function") 

(Nomenclature "High Pressure Oxidizer Turbine Turbine Torque Characteristics") 
(Input-parameter ) ) 

(Deframe Trot2 

(Class TURBINE) 
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(Unit "Table-Function") 

(Nomenclature "Low Pressure Oxidizer Turbine Turbine Torque Characteristics”) 
(Input-parameter ) ) 

(Deframe Tpr 

(Class TURBINE) 

(Unit "Table-Function") 

(Nomenclature "Pressure Ratio") 

(Input-parameter ) ) 


ORIGINAL PAGE 13 
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ORIGINAL PAGE IQ 
OF POOR QUALITY 


(Deframe Sfl 

(Class TURBOPUMP) 

(Unit "Rotational-Speed") 

(Nomenclature "Low Pressure Fuel Turbopump Rotational Speed") 
(Input-parameter B14 TQftl-TQfpl SflO) ) 

(Deframe Sf2 

(Class TURBOPUMP) 

(Unit "Rotational-Speed") 

(Nomenclature "High Pressure Fuel Turbopump Rotational Speed") 
(Input -parameter B23 TQft2-TQfp2 Sf20) ) 

(Deframe Sol 

(Class TURBOPUMP) 

(Unit "Rotational-Speed") 

(Nomenclature "Low Pressure Oxidizer Turbopump Rotational Speed") 
(Input-parameter B35 TQotl-TQopl SolO)) 

(Deframe So2 

(Class TURBOPUMP) 

(Unit "Rotational-Speed") 

(Nomenclature "High Pressure Oxidizer Turbopump Rotational Speed") 
(Input-parameter B42 TQot2-TQop2-TQop3 So20) ) 

(Deframe SflO 

(Class TURBOPUMP) 

(Unit "Initial-Condition") 

(Nomenclature "Sfl when t=0") 

(Input-parameter 1638.0)) 

(Deframe Sf20 

(Class TURBOPUMP) 

(Unit "Initial-Condition") 

(Nomenclature "Sf2 when t=0") 

(Input-parameter 3922.0)) 

(Deframe SolO 

(Class TURBOPUMP) 

(Unit "Initial-Condition") 

(Nomenclature "Sol when t=0") 

(Input-parameter 570.9)) 

(Deframe So20 

(Class TURBOPUMP) 

(Unit "Initial-Condition") 

(Nomenclature ”So2 when t=0") 

(Input-parameter 3262.0)) 
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(De frame Amfv 

(Class VALVE) 

(Unit "Area") 

(Nomenclature "Main Fuel Valve Area") 

(Input-parameter TQmfv Xmfv) ) 

(Deframe Amov 

(Class VALVE) 

(Unit "Area") 

(Nomenclature "Main Oxidizer Valve Area") 
(Input-parameter TQmov Xmov) ) 

(Deframe Afpv 

(Class VALVE) 

(Unit "Area") 

(Nomenclature "Fuel Preburner Oxidizer Valve Area") 
(Input-parameter TQfpv Xfpv) ) 

(Deframe Aopv 

(Class VALVE) 

(Unit "Area") 

(Nomenclature "Oxidizer Preburner Oxidizer Valve Area") 
(Input-parameter TQopv Xopv) ) 


ORIGINAL PAGE IS 
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System Generated Equations 


Al:>PSAM>equation>psam-equation.eq.36 


10/03/8815:28:59 Page 1 


Component Processor Element Pump 


Component Processor Element Turbine 


S02 = B42 T (TQOT2 - TQOP3 - TQOP2 ) + SO20 

SOI = B35 T (TQOT1 - TQOP1) + SO10 

SF2 = B23 T (TQFT2 - TQFP2) + SF20 

SF1 = B14 T (TQFT1 - TQFP1) + SF10 


B48 POP TPR 

DWOT2 = SQRT ( ) 

TOP 


POD2 - POD1 

DWOT1 = SQRT ( ) 

ROT1 + B37 

DWFT2 = SQRT (B24 PFP TPR) 

DWFT1 = SQRT (B1 7 (P5 - PFI) R5) 

B47 S02 

UOT2 = 

SQRT (ROP TOP) 

B38 SOI 

FOTl = 

SQRT (DWOT1) 

B25 SF2 

UFT2 = 

SQRT <TFP) 

B15 SF1 

UFTl - 

SQRT (T5 ) 

TQOT2 = B29 POP TROT2 
2 

TQOT1 = B36 DWOT1 TROT1 
TQFT2 = B26 PFP TRFT2 
TQFT1 = B16 P5 TRFT1 

Component Processor Element Turbopump 

/ 

[ 1 

S02 (T) = B42 (I dT> {TROT2 - TROP3 - TROP2) + S0O2 (T) 

] G(T) 

/ 

/ 

C 1 

SOI (T) = B35 {I dT) (TROT1 - TROP1) + SOOl (T) 

ORIGINAL PAGE IS > G(T > 

OF POOR QUALITY 

/ 

[ 1 

SF2 (T) = B23 (I dT) (TRFT2 - TRFP2 ) + S0F2 (T) 

] G<T) 

/ 

/ 

t 1 

SF1 (T) = B14 (I dT) (TRFT1 - TRFP1 ) + S0F1 (T) 

] G(T) 

/ 
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A1 : >PS AM>equation>psam-equation.eq.36 


10/03/8815:28:59 Page 2 


Component Processor Element Valve 


Component Processor Element Flow 


00 (T) 


00 (T) 


AOPV = TQOPV XOPV 
AFPV = TQFPV XFPV 
AMOV = TQMOV XMOV 
AMFV = TQMFV XMFV 

DWOPOI = DWOPO EOPO 

DWFPOI = DWFPO EFPO 

2 2 
AB B61 DWOPO (T) 


ORIGINAL PAG.E fS 
OF POOR QUALITY 


DWOPO (T) = B59 I (- B60 DWOPOI (T) - 


+ PPOS - POP) dT + DWOP 


/ 

[ 2 
DWFPO (T) - B52 I {- B54 DWFPOI (T) - 


2 2 
AB B53 DWFPO (T) 


+ PPOS - PFP) dT + DWFP 


/ 

B57 [ 2 

DWOPF(T) = B <- - POP + P9) I DWOPF (T) dT + DWOO (T) 

R9 ] 

/ 

/ 

[ 2 


B50 


DWFPF(T) = B (- - PFP + P9) I DWFPF (T) dT + DWFO (T) 

R9 ] 

/ 


Component Processor Element Combuster 


TOP = TTFOP + B113 T9 
TFP = TTPFP + B112 T9 
DWOPOI 

DWOPOI + DW1 
DWFPOI 


FOP = 


FFP 


DWFPOI + DW1 


/ 

[ 

POP = B63 I (- B107 DWOT2 (T) + DWOPOI (T) + DWOPF (T) ) dT + POPO 

] 

/ 

/ 

c 

PFP = B56 I (- Bill DWFT2 (T) + DWFPOI (T) + DWFPF (T) ) dT + PFPO 

] 

/ 

Component Processor Element Cooling 

/ 

[ 

TW24 = TW240 - B93 I DQW24 (T) dT 

] 

/ 
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A1 : >PS AM>equation>psam-equation.eq.36 


10/03/8815:28:59 Page 3 


/ 

[ 

TW1 4 = TW140 + B92 I (DQTC4 (T) - DQW14(T)) dT 

] 

/ 

/ 

[ 

ORIGINAL PAGE IS TW25 = TW250 - B83 I DQW25 (T) dT 

OF POOR QUALITY \ 

/ 

[ 

TW15 = TW150 + B82 I (DQTC5 (T) - DQW15(T)) dT 

] 

/ 

0.8 

DQTC4 = B91 DWCN (TC - TW14) 

0.8 

DQTC5 = B81 DWCN (TC - TW15) 

0.8 

DQW24 = B90 DWFN (0.002 T4 + 1.0) (TW24 - T4) 

0.8 

DQW14 = B89 DWFN (0.002 T4 + 1.0) (TW14 - T4) 

0.8 

DQW25 = B80 DWMC (0.002 T5 + 1.0) (TW25 - T5) 

0.8 

DQW15 = B79 DWMC (0.002 T5 + 1.0) (TW15 - T5) 

B88 P 4 

T4 

R4 

B78 P5 

T5 = 

R5 

/ 

[ 

R4 = B87 I (DWFN (T) - DW4 (T) ) dT + R40 

3 

/ 

/ 

[ 

R5 = B77 I (DWMC(T) - DW5 (T) ) dT + R50 
3 
/ 

/ 

[ 

DW4 (T) = B94 I { 

3 

/ 

2 

B76 DWMC (T) 

+ PMFVD - P 5 ) dT + DWMCO 

R5 


H4 = B102 T 4 
H5 = B101 T5 

/ 

[ 

P 4 = B84 I (- H3 DWFN (T) + H4 DW4 (T) + DQW2 4 (T) + DQW14 (T) ) dT + P40 

I 


/ 

[ 

DWMC (T ) = B75 I (- 
3 

/ 


2 

B95 DW4 (T) 

-P9+P4) dT + DW40 

R4 
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Al:>PSAM>equation>psam-equation.eq.36 


10/03/88 15:28:59 Page 4 


ORIGINAL IS ' 

OF POOR QUALITY { 

P5 = I (H3 DWMC(T) - H5 DW5 (T) + DQW25 (T) + DQW15 (T) ) dT + P50 
3 
/ 

Component Processor Element flowmeter 

TOD2 = TOS + 0.00547 S02 
TFP1 = TFS + 0.0015789 SF1 


OSS - 164.0) ) ] 

1 


POFM = POD2 


2 

(DWOTPR + DWOP3 + DWMOV) RFOD2 
ROD 2 
2 


DWFD2 RFFCD 

PFFM = PFD1 - 

RFD1 

QOFM = [(0.0209912 S02 + 6.198 91 (DWOTPR + DWOP3 + DWMOV) + 17.5339) (1.0 


0.002789 (T 


QF 


[(0.20206 SF1 + 101.323 DWFD2 - 204.852) (1.0 


0.008783 (TFS - 37.0) ) 


Component Processor Element Preburner Supply Duct 

DWFNBP + DW4 

R9 * [ ] 

DW4 

+ BIOS DWFNBP 

R4 

B100 P 9 

T9 

R9 

/ 

[ 

PPOS = B49 I (- DWOPO(T) + DWOP3 (T) - DWFPO(T) + DW(T)) dT + PPOSO 

3 

/ 

/ 

C 

P 9 = B9 6 I (- DWOPF(T) - DWFPF (T) + 2 DWFNBP (T) ) dT + P90 

3 

/ 

Component Processor Element High Pressure Oxidizer Pump Preburner Boost-stage 

/ 

[ 2 

DWOP3 (T) = B1 I (- B2 DWOP3 (T) + P3 - P) dT + DW0 
3 
/ 

2 

POD3 = B F S TT + P2 
2 

TOP 3 = B44 F S TT 

B DW 

FOP 3 = 

S 

Component Processor Element Main Chamber Injector, Fuel, Hot Gas and Oxidizer 


EOT 


TEOI WOI 
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A1 : >PS AM>equation>psam-equation.eq.36 


10/03/8815:28:59 Page 5 


ORIGINAL PAGE IS 
OF POOR QUALITY 


/ 

c 

WOI = WOI + I (DWMOV(T) - DWOI (T) ) dT 

] 

/ 


TFI 


DWOI = DWMOV EOI 
B67 TOP + B66 TFP + B68 T5 
B64 PC PFI TPR 


DWFI 


0.5 

PPI TFI 


/ 

[ 

PFI = B 64 I (DWOT2 (T) + DWFT2 (T) + DWFT1 (T) - B108 DWFI(T)) dT + 

] 

/ 

Component Processor Element Main Thrust Chamber 

/ 

t 

PC = B72 I (DWOI(T) + B109 DWFI (T) - DWCN (T) ) dT + PCO 

] 

/ 

PC = B114 PC 

Component Processor Element Nozzle Bypass Element 

DWCCV = DWFNBP 


/ 

[ B99 DWFNBP (T) 

DWFNBP = B98 I (- + PMFVD - P9) dT + DWFNBP 0 

] 2 
/ ACCV 


PFIO 
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System Generated FORTRAN Code 



A1 : >PS AM>FORTRAN>psatn-Code.for.21 


10/03/8815:36:34 Page 1 


Component Processor Element Pump 

TORQUE = B41*S02**2*TR0P2 
ROP1 - B34*S01**2*TR0P1 
RFP2 = B20*SF2**2*TRFP2 
RFP1 = B13*SF1**2*TRFP1 
POD2 = B40*SO2**2*TPOP2+PODl 
POD1 = B28*S01**2*TPOP1+POS 
PFD2 = B19*SF2**2*TPFP2+PFD1 
PFD1 = B12*SF1**2*TPFP1+PFS 
FOP2 - B39* (DWOT1+DWOP3+DWMOV) /S02 
FDP1 = B27* (DWOP3+DWMOV) /SOI 
FFP2 - B18*DWFD2/SF2 
FFP1 - B11*DWFD1/SF1 
Component Processor Element Turbine 


ORIGINAL PAGE IS 
OF POOR QUALITY 


DWOT2 = SQRT (B48*POP*TPR/TOP) 

DWOT1 - SQRT ( (POD2-POD1) / (ROT1+B37) ) 

DWFT2 = SQRT (B24*PFP*TPR) 

DWFT1 - SQRT (B17* (P5-PFI) *R5) 

UOT2 = B47*S02/ SQRT (ROP*TOP) 

FOT1 = B38*S01/ SQRT (DWOT1 ) 

UFT2 = B25*SF2/SQRT (TFP) 

UFT1 - B15*SF1/ SQRT (T5) 

TQOT2 = B2 9* POP* TROT 2 
TQOT1 = B36*DW0T1**2*TR0T1 
TQFT2 = B2 6*PFP*TRFT2 
TQFT1 - B16*P5*TRFT1 
Component Processor Element Turbopump 

S02 = B42*T* (TQOT2-TQOP3-TQOP2 ) +S020 
SOI - B35*T* (TQOT1-TQOP1) +SO10 
SF2 - B23*T* (TQFT2-TQFP2) +SF20 
SF1 - B14*T* (TQFT1-TQFP1 ) +SF10 
Component Processor Element Valve 

AOPV = TQOPV*XOPV 
AFPV - TQFPV*XFPV 
AMOV = TQMOV*XMOV 
AMFV - TQMFV*XMFV 
Component Processor Element Flow 

DWOPOI - DWOPO*EOPO 
DWFPOI - DWFPO*EFPO 

DWOPO (T) = B59* ' INTEGRATE (-B60*DWOPOI (T) **2-AB**2*B61*DWOPO (T) **2/ 
1 A**2+PPOS-POP, T) +DWOPOO (T) 

DWFPO (T) - B52* ' INTEGRATE (-B54*DWFPOI (T) **2-AB**2*B53*DWFPO (T) **2/ 
1 A**2+PPOS-PFP, T) +DWFPOO (T) 

DWOPF(T) - B* (-B57/R9-POP+P9) *' INTEGRATE (DWOPF (T) **2, T) +DWOO (T) 
DWFPF(T) = B* (-B50/R9-PFP+P9) *' INTEGRATE (DWFPF (T) **2, T) +DWF0 (T) 
Component Processor Element Combuster 

TOP = TTFOP+B113*T9 
TFP = TTPFP+B112*T9 
FOP - DWOPOI/ (DWOPOI+DW1) 

FFP = DWFPOI/ (DWFPOI+DW1) 

POP = B63* ' INTEGRATE (-B107*DWOT2 (T) +DWOPOI (T) +DWOPF (T) , T) +POP0 
PFP = B56* ' INTEGRATE (-B111*DWFT2 (T) +DWFPOI (T) +DWFPF (T) , T) +PFP0 
Component Processor Element Cooling 
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TW24 = TW240-B93*' INTEGRATE (DQW24 (T) , T) 


Al:>PSAM>FORTRAN>psam-Code.for.21 


10/03/8815:36:34 Page 2 


TW14 - TW140+B92*' INTEGRATE (DQTC4 (T)-DQW14 (T) , T) 

TW25 - TW250-B83*' INTEGRATE (DQW25 (T) , T) 

TW15 - TW150+B82*' INTEGRATE (DQTC5 (T)-DQW15 (T) ,T) ORIGINAL 

DQTC4 = B91*DWCN**0 . 8* (TC-TW14) _ L IS 

DQTC5 = B81*DWCN**0 . 8* (TC-TW15) ut H00R Quality 

DQW24 - B90*DWFN**0 .8* (0.002*T4 + 1.0)* (TW24-T4 ) 

DQW14 - B8 9*DWFN**0 ,8*(0.002*T4+1.0)* (TW14-T4 ) 

DQW25 = B80*DWMC**0 . 8* (0 . 002*T5+1 . 0) * (TW25-T5) 

DQW15 - B7 9*DWMC**0 .8*(0.002*T5 + 1.0)* (TW15-T5 ) 

T4 - B88*P4/R4 
T5 - B78*P5/R5 

R4 = B87*' INTEGRATE (DWFN (T) -DW4 (T) , T) +R40 

R5 - B77*' INTEGRATE (DWMC (T) -DW5 (T) , T) +R50 

DW4 (T) = B94*' INTEGRATE (-B95*DW4 (T) **2/R4-P9+P4 , T) +DW40 

DWMC(T) - B75* ' INTEGRATE (-B7 6*DWMC (T) **2 /R5+PMFVD-P5, T) +DWMC0 

H4 - B102*T4 

H5 = B101*T5 

P4 - B84* ' INTEGRATE (-H3*DWFN (T ) +H4*DW4 (T) +DQW24 (T) +DQW14 (T) , T) +P40 
P5 - ' INTEGRATE (H3*DWMC (T) -H5*DW5 (T) +DQW25 (T) +DQW15 (T) , T) +P50 
Component Processor Element flowmeter 

TOD2 - TOS+O .00547*S02 
TFP1 = TFS+0 . 0015789*SF1 

POFM = POD2- (DWOTPR+DWOP3+DWMOV) **2 *RFOD2/ROD2 
PFFM - PFD1-DWFD2**2 *RFFCD/RFD1 

QOFM = [ (0 . 0209912*SO2 + 6. 19891* (DWOTPR+DWOP3+DWMOV) +17.5339) * (1.0- 
1 0.002789* (TOSS-164.0) ) ] 

QF = [ (0 . 20206*SF1+101 . 323*DWFD2-204 . 852 ) * (1.0-0.008783* (TFS-37.0) 

1 ) ] 

Component Processor Element Preburner Supply Duct 

R9 - [ (DWFNBP+DW4) / (DW4/R4+B105*DWFNBP) ] 

T9 = B100*P9/R9 

PPOS = B49*' INTEGRATE (-DWOPO (T) +DWOP3 (T) -DWFPO (T) +DW (T) , T) +PPOS0 
P9 = B96* ' INTEGRATE (-DWOPF (T) -DWFPF (T) +2*DWFNBP (T) , T) +P90 
Component Processor Element High Pressure Oxidizer Pump Preburner Boost-stage 

DWOP3 (T) - Bl*' INTEGRATE (-B2*DWOP3 (T) **2+P3-P, T) +DW0 
POD3 = B*F*S**2*TT+P2 
TOP 3 - B44*F*S**2*TT 
FOP 3 - B*DW/S 

Component Processor Element Main Chamber Injector, Fuel, Hot Gas and Oxidizer 
EOI - TEOI*WOI 

WOI - WOI+' INTEGRATE (DWMOV(T) -DWOI (T) , T) 

DWOI - DWMOV*EOI 

TFI - B67*TOP+B66*TFP+B68*T5 

DWFI - B64*PC*PFI*TPR/ (PPI*TFI**0 . 5) 

PFI - B64*' INTEGRATE (DWOT2 (T)+DWFT2 (T) +DWFT1 (T) -B108*DWFI (T) ,T)+PF 
1 10 

Component Processor Element Main Thrust Chamber 

PC - B72*' INTEGRATE (DWOI (T) +B109*DWFI (T) -DWCN (T) , T) +PC0 
PC - B114*PC 

Component Processor Element Nozzle Bypass Element 
DWCCV - DWFNBP 

DWFNBP - B98*' INTEGRATE (-B99*DWFNBP (T) /ACCV**2+PMFVD-P9, T) +DWFNBP0 
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Documentation for Component Processure PUMP 

1. Low Pressure Fuel Pump 

2. High Pressure Fuel Pump 

3. Low Pressure Oxidizer Pump 

4. High Pressure Oxidizer Pump 
Generic Equations 

F ■■ Bi * [ DW / S] ; Flow variable 

P = Ps + Bj * (S)-2 * Tp (F ) ; Total pressure 

R = Bk * (S) A 2 * Tr(F) ; Torque 

EOPOO : Eopo when t=0 

EFPOO : EFpo when t=0 

DWOPOO : DWopo when t=0 

DWFPOO : DWfpo when t=0 

DWOPFO : DWopf when t=0 

DWFPFO : DWfpf when t-0 

TROP2 : High Pressure Oxidizer Pump Pump Torque Characteristics 

TORQUE : High Pressure Oxidizer Pump "Torque" 

TROPl : Low Pressure Oxidizer Pump Pump Torque Characteristics 

TRFP2 : High Pressure Fuel Pump Pump Pump Torque Characteristics 

TRFP1 : Low Pressure Fuel Pump Pump Torque Characteristics 

TP0P2 : High Pressure Oxidizer Pump Pump Pressure Rise Characteristics 

TPOPl : Low Pressure Oxidizer Pump 

TPFP2 : High Pressure Fuel Pump Pump Pressure Rise Characteristics 
TPFPl : Low Pressure Fuel Pump Pump Pressure Rise Characteristics 
ROPl : Low Pressure Oxidizer Pump 
RFP2 : High Pressure Fuel Pump 
RFPl : Low Pressure Fuel Pump 

POD2 : High Pressure Oxidizer Pump Total Pressure 

PODl : Low Pressure Oxidizer Pump Total Pressure 

PFD2 : High Pressure Fuel Pump Total Pressure 

PFD1 : Low Pressure Fuel Pump Discharge Total pressure 

FOP2 : High Pressure Oxidizer Pump Flow variable 

FDPl : Low Pressure Oxidizer Pump Flow variable 

FFP2 : High Pressure Fuel Pump Flow Variable 


, in^3 
, Psia 
, in-lb 



A1 : >PS AM>doc>psam-doc.doc.38 


10/03/88 14:08:46 Page 2 


FFPl : Low Pressure Fuel Pump Flow variable 
*******★*■*■****★★•*******★*■*■*■■*★*★*★**■■*■■*■★★**★★★★*■*★*★ 


Documentation for TURBINE section 

1. Low Pressure Fuel Turbine 

2. High Pressure Fuel Turbine 

3. Low Pressure Oxidizer Turbine 

4. High Pressure Oxidizer Turbine 

Generic Equations 

R = Bi * P * Tr ; Torque , in-lb f 

Ts = Bj * S / (T A ( 1/2 ) ) ; Speed Parameter 

R = Bk * (S) A 2 * Tr (F) ; Weight Flowrate , lb m/s ec 

TPR : Pressure Ratio 

TROT 2 : Low Pressure Oxidizer Turbine Turbine Torque Characteristics 

TROT1 : High Pressure Oxidizer Turbine Turbine Torque Characteristics 

TRFT2 : High Pressure Fuel Turbine Turbine Torque Characteristics 

TRFTl : Low Pressure Fuel Turbine Turbine Torque Characteristics 

DWOT2 : High Pressure Oxidizer Turbine Weight Flowrate 

DWOT1 : Low Pressure Oxidizer Turbine Weight Flowrate 

DWFT2 : High Pressure Fuel Turbine Weight Flowrate 

DWFTl : Low Pressure Fuel Turbine Weight Flowrate 

UOT2 : High Pressure Oxidizer Turbine Speed Parameter 

FOTl : Low Pressure Oxidizer Turbine Flow Variable 

UFT2 : High Pressure Fuel Turbine Speed Parameter 

UFT1 : Low Pressure Fuel Turbine Speed Parameter 

TQOT2 : High Pressure Oxidizer Turbine 

TQOT1 : Low Pressure Oxidizer Turbine Torque 

TQFT2 : High Pressure Fuel Turbine Torque 

TQFT1 : Low Pressure Fuel Turbine Torque 

Documentationf or TURBOPUMP section 

1. Low Pressure Fuel Turbopump 

2. High Pressure Fuel Turbopump 

3. Low Pressure Oxidizer Turbopump nDrfMMJll oa 

4. High Pressure Oxidizer Turbopump UKlGIINIAL PAGE IS 

Generic Equations OF POOR QUALITY 

S = (l/G) * 'Integrate (Tq*temp (t) , t) +S0 ; Speed 
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SO20 : So2 when t-0 
SO10 : Sol when t=»0 
SF20 : Sf 2 when t-0 
SF10 : Sfl when t-0 

S02 : High Pressure Oxidizer Turbopump Rotational Speed 
SOI : Low Pressure Oxidizer Turbopump Rotational Speed 
SF2 : High Pressure Fuel Turbopump Rotational Speed 
SF1 : Low Pressure Fuel Turbopump Rotational Speed 
************************************************** 
************************************************** 


Create Knowledge Base for VALVE section 

1 . Main Fuel Valve 

2. Main Oxidizer Valve 

3. Fuel Preburner Oxidizer Valve 

4. Oxidizer Preburner Oxidizer Valve 
Generic Equations 

A “ Ab * (T * Th | Thb) ) 

Input-parameter : Ab T Th Thb 

Ab: Average of area, T * Th | Thb =Valve Effective Area 
Pit *= Pi - RF *DW A 2 | RHO 

Input-parameter ; Pi RF DW RHO 

Pit: Inlet Total Pressure, Pi: Inlet Pressure 
RF: Resistance to Flow, DW: Weight Flowrate, Rho: Density 

Pis - Pit - DW A 2 | [(772.8)* (Ai A 2) * RHO] 

Input-parameter : Pit DW Ai RHO 
dP = B * (RHOb | RHO) (DW | (AjAb) A 2 - RF * DW | RHO 
Input-parameter : B Rhob Rho DW A Ab RF 
AOPV : Oxidizer Preburner Oxidizer Valve Area 


AFPV : Fuel Preburner Oxidizer Valve Area 
AMOV : Main Oxidizer Valve Area 
AMFV : Main Fuel Valve Area 
************************************************** 




Create Knowledge Base for FLOW section 

1. Fuel Preburner Fuel Flow 

2. Fuel Prebumer Oxidizer Flow 

3. Oxidizer Preburner Fuel Flow 

4. Oxidizer Preburner Oxidizer Flow 
Generic Equations 

DWf - Bi INT { P { 9) -P- [B j /R ( 9 ) ] (DWf) A 2}dt +[DWf]t=0 ;Fuel Flow 

DWo = Bk INTfPpos - P - Be[DWo/ <A/A' ) ] A 2 - Bm(DWi) A 2}dt + [DWo] t-0 /Oxidizer Flow 
Eo == Bn INT {DWo - DWi) }dt +[Eo]t=0 ; Variable in injector Priming function 
DWi = DWo [Eo] 

DWOPOI : Oxidizer Preburner Oxidizer Injector Weight Flowrate 
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DWFPOI : Fuel Preburner Oxidizer Injector Weight Flowrate 
EOPO : Oxidizer Preburner Oxidizer Variable in injector Priming 
EFPO : Fuel Preburner Oxidizer Variable in injhector Priming 
DWOPO : Oxidizer Preburner Oxidizer Weight Flowrate 
DWFPO : Fuel Preburner Oxidizer Weight Flowrate 
DWOPF : Oxidizer Preburner Fuel Flow 
DWFPF : Fuel Preburner Fuel Flow 
************************************************** 
****•*★★★**★*★**★*★********★*****★***★*★★********** 

Create Knowledge Base for COMBUSTER section 

1. Fuel Preburner COMBUSTER 

2. Oxidizer Preburner COMBUSTER 
Generic Equations 

P = Bi * Integrate { [DW1 + DW2 - Bj * DW3],t} + P0 ; Total Pressure, Psia 
F - DW2/<DW1+ DW2) /Fraction 
T = Tqt (F) + Bk * T9 /Temperature 

TTFOP : Oxidizer Preburner Hot Gas Temperature Function 

TTFFP : Fuel Preburner Hot Gass Temperature Function 

POPO : Pop when t=0 

PFPO : Pfp when t=0 

TOP : Oxidizer Preburner Combuster Temperature 
TFP : Fuel Preburner Combuster Temperature 
FOP : Oxidizer Preburner Flowrate Fraction 
FFP : Fuel Preburner Flowrate Fraction 
POP : Oxidizer Preburner Total Pressure 
PFP : Fuel Preburner Total Pressure 


ORIGINAL page is 

OF POOR QUALITY 


Create Knowledge Base for COOLING section 

1 . Main Chamber Regenerative Cooling Element 

2. Nozzle Regenerative Cooling Element 
Generic Equations 

P * B * integrate { [DQwl + DQw2 + H3 * DWi - H * DW],t} + P0 / Total Pressure 
H - b * T / Specific Enthalpy 

DW «Bi * integrate { [P - P5 - (B j/R5) * DW A 2],t} +DW0 / Weight Flowrate 
R ** B * integrate { [DWi - DW],t} + R0 ; Density 
T5 - B * (P / R) / Temperature 

DQwl - b *[1.0 + 0.002 * T5] * [Twl5 - T5]* DWmc A 0.8 /Heat transfer rate 
DQw2 «= B * [1.0 + 0.002 * T5]*[Tw25 - T5]*DWmc A 8 / Heat Transfer Rate 
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DQtc - B * [Tc - Twl5 ] * [DWcn] A 0 . 8 

Twl « B * integrate {[DQtc - DQwl],t} + TwlO ; Hot Gas Wall Temperature 
Tw2 “ B * integrate {[-DQw2] / t} + Tw20 ; Ambient thrust chamber Temperature 
TW240 : Tw2 4 when t-0 

TW140 : Twl 4 when t-0 

TW250 : Tw2 5 when t=0 

TW150 : TWl 5 when t-0 

R40 : R4 when t=0 

R50 : R5 when t-0 

DW40 : DW4 when t-0 

DWMCO : DWmc when t-0 

P40 : P4 when t-0 

P 50 : P5 when t-0 

TW24 : Ambient Thrust Temperature 

TW1 4 : Hot Gas temperature 

TW25 : Ambient thrust chamber Temperature 

TW15 : Hot Gas Temperature 

DQTC4 : Nozzle Regenerative Heat Transfer Rate TC 

DQTC5 : Main Chamber Cooling Element Heat Transfer Rate TC 

DQW24 : Nozzle Regenrative Ambient Thrust Heat Transfer Rate 

DQW14 : Nozzle Regenerative Hot Gas Heat Transfer Rate 

DQW25 : Main Chamber Ambient Thrust Chamber Heat Transfer Rate 

DQW15 : Main Chamber Hot Gas Heat Transfer Rate 

T4 : Nozzle Regenerative Cooling Temperature 

T5 : Main Chamber Heat Exchanger Temperature 

R4 : Nozzle Regenerative Element Density 

R5 : Main Chamber Element Density 

DW4 : Nozzle Regenerative Cooling Weight Flowrate 
DW5 : Main Chamber Heat Weight Flowrate 

H4 : Nozzle Regenerative Cooling Element Specific Enthalpy 

H5 : Main Chamber Heat Exchanger Specific Enthalpy 

P4 : Nozzle Regenerative Exchanger Element Total Pressure 


P5 : Main Chamber Heat Exchanger Total Pressure 
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************************************************** 

************************************************** 

Create Knowledge Base for FLOWMETER section 

1. Fuel Flowmeter 

2. Oxidizer Flowmeter 
Generic Equations 

Q = P -(RF * DW A 2)/Rb ; Volumetric Flowrate , GPM 
P = Pf - (RF * DW A 2)/R ; Total Pressure , Psia 
T = Ti +B * S ; Temperature 

TOD2 : Low pressure Oxidizer pump Flowmeter Temperature 
TFPl : Low Pressure Fuel Pump Flowmeter Temperature 
PFFM : Fuel Flowmeter Total Pressure 
QF : Fuel Volumetric Flowrate 
************************************************** 
************************************************** 

Create Knowledge Base for PREBURNER- SUPPLY-DUCT section 

1 . Preburner Fuel Supply Duct 

2. Preburner Oxidizer Supply Duct 
Generic Equations 

P = B * integrate! [DW] , t} + P0 ; Total Pressure , Psia 
T = B * [P/R] ; Temperature 

R * [DW1 + DW2] / [ (DW2/R) + B * DW3 ] ; Density 

PPOS0 : Ppos when t=0 

P 90 : P9 when t=0 

R9 : Preburner Fuel Supply Element Density 
T9 : Preburner Fuel Supply Element Temperature 
PPOS : Preburner Oxidizer Supply Total Pressure 
P 9 : Preburner Fuel Supply Element Total Pressure 
************************************************** 


************************************************** 

Create Knowledge Base for BOOST-STAGE section 

1. High Pressure Oxidizer boostpump Flow variable 

2. High Pressure Oxidizer boostpump rotational Torque 

3. High Pressure Oxidizer Boostpump pressure 

4. High Pressure Oxidizer Boostpump weight flowrate 
Generic Equations 

F - B * (DW / S)) ; Flow variable , in A 3 

T = B * (S A 2 ) * TT * (F)) ; Torque , in-lb 

P3 = B * (S A 2 ) * TT * F + P2) ; Pressure 

DW3 - B1 * integrate [ (P3 - P - B2 A 2) , t] + DWO); Weight Flowrate, lbm/sec 

DWOP30 : DWop3 when t=0 

TP0P3 : high Pressure Oxidizer boostpump Pressure Rise Character 
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TROP3 : High Pressure Oxidizer Boostpump Torque 
DWOP3 : High Pressure Oxidizer Boostpump Weight Flowrate 
POD3 : High Pressure Oxidizer Boostpump Pressure 
TOP3 : High Pressure Oxidizer Boostpump Torque 
FOP3 : High Pressure Oxidizer boostpump Flow variable 
************************************************** 
************************************************** 

Create Knowledge Base for BOOST-STAGE section 

1. High Pressure Oxidizer boostpump Flow variable 

2. High Pressure Oxidizer boostpump rotational Torque 

3. High Pressure Oxidizer Boostpump pressure 

4. High Pressure Oxidizer Boostpump weight flowrate 
Generic Equations 

F “ B * (DW / S) ) ; Flow variable , in A 3 

T = B * (S A 2 ) * TT * (F)) ; Torque , in-lb 

P3 * B * (S A 2 ) * TT * F + P2 ) ; Pressure 

DW3 = B1 * integrate [ (P3 - P - B2 A 2) , t] + DWO); Weight Flowrate, lbm/sec 

PFIO : Pfi when t=0 

EOI : Main Chamber injector. Oxidizer variable in injector priming function 

WOI : Main Chamber Injector, Oxidizer Weight 

DWOI : Main Chamber injector. Oxidizer Weight Flowrate 

TFI : Main Chamber Temperature 

DWFI : Main Chamber injector weight flowrate 

PFI : Main Chamber Injector, fuel and Hot Gas Total Pressure 
************************************************** 
************************************************** 


Create Knowledge Base for MAIN- THRUST-CHAMBER section 

1. Combuster Injector end static Total Pressure 

2. Main thrust Chamber Total Pressure 

Generic Equations 

P = b * P ; Combuster Injection Total pressure , Psia 

P - B1 * 'integrate [ (B2 * DWl + DW2 - DW3),t] + P0) ; Main Thrust Chamber Total Pr 
essure, Psia 
PCO : Pc when t=0 

PC : Main thrust chamber Total Pressure 
PCIES : Combuster Injector end static Total Pressure 
**********************************^*************** 
★★★a********************************************** 


Create Knowledge Base for NOZZLE-BYPASS section 
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1. Nozzle Heat exchanger bypass Weight Flowrate 

2 . Coolant Control valve 
Generic Equations 

DW ~ Bl * ' integrate ( [P - P9 - B2 * ( DW1/ (A/Ab) A 2 ] , t ) + DWO ; Nozzle Heate weight 

flowrate , lbm/sec 

DW = DWl ; Coolant Control Valve , lbm/sec 
DWFNBPO : DWfnbp when t=0 

DWCCV : Coolant Control Valve 

DWFNBP : Nozzle Heat Exchanger bypass weight flowrate 





Example for the Robot Manipulators 

Example steps to solve the end position. 

1 . Define rules. sin 2 0 t = 1 — cos 2 0j 

SSj = sin 8 j i is from 1 to 6 

CCj = cos 0j i is from 1 to 6 

2. Define matrices A|, A 2 , A 3 , A 4 , A 5 , A 6 , T 6 . 

3. T 5 = A 1 A 2 A 3 A 4 A 5 A^ 

T 6 is the end position and A* is the matrix for each joint. 

T 6 is the product of A| ... A 6 

4. A| _l Tg = A 2 A 3 A 4 Ag A^ 

A 1 1 A 2 _i T 6 = A 3 A 4 A 5 A 6 

A 1 1 A 2 1 ... A 5 1 T 5 = A 6 

5. Compare the above equations and solve 
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End position T 6 


n x o x a x 

n y o y a y 

n* o 7 a z 

0 0 0 


6. The end position is a product of the A matrices 

T 6 = A | A2 • A3 • A4 • A 5 • A 6 

7 . Matrix inverse, matrix multiplication 


T,6 

= 

A,' 1 

t 6 

= a 2 

> 

> 

• As 

■ a 6 

T 2 6 

= 

A,' 1 

a 2 - 

1 T 6 = 

- A3 • A4 

• As ' 

■ a 6 

T36 

= 

A,' 1 

a 2 - 

1 A3- 1 

t 6 = a 4 

• As • 

a 6 

T46 

= 

A,' 1 

A 2 - 

1 A3- 1 

a 4 1 t 6 - 

: As • 

Aft 

T 5 6 


A,' 1 

a 2 - 

1 A3- 1 

A4 1 A 5 1 

t 6 = 

A 6 


8. Evaluating equations to obtain the end position 
n x , n y , n 7 , o x , o y , o 7 , a x , a y , a 7 , p x , p y , p 7 


ORIGINAL PAGE IS 

OF POOR QUALITY 
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(C28) sa3i: rat simp (a3i) ; 
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